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Motivated by several recent data, we test the QCD spectral sum rules (QSSR) predictions based on different proposals ( qq , 
qqqq, and gluonium) for the nature of scalar mesons. In the 7 = 1 and 1/2 channels, the unusual wrong splitting between 
the ao(980) and k(900) and the ao(980) width can be understood from QSSR within a qq assignement. However, none of 
the qq and qqqq results can explain the large k width, which may suggest that it can result from a strong interference with 
non-resonant backgrounds. In the 1 = 0 channel, QSSR and some low-energy theorems (LET) require the existence of a 
low mass gluonium <tb( 1 GeV) coupled strongly to Goldstone boson pairs which plays in the U(l)v channel, a similar role 
than the r\ for the value of the 1/(1) a topological charge. The observed cr(600) and /o(980) mesons result from a maximal 
mixing between the gluonium ob and qq( 1 GeV) mesons, a mixing scheme which passes several experimental tests. OZI 
violating J/%j) —> (f)n + n~, D s —> 3n decays and J/ip — > jS glueball filter processes may indicate that the /o(1500), /o(1710) 
and /o(1790) have significant gluonium component in their wave functions, while the /o(1370) is mostly qq. Tests of these 
results can be provided by the measurements of the pure gluonium r/g and 4-7T specific U(1)a decay channels. 


1. Introduction 

The nature of scalar mesons continues to be an intrigu¬ 
ing problem in QCD. Experimentally, there are well es¬ 
tablished scalar mesons with isospin 7=1, the cio(980), 
ao(1450) with isospin I = 1/2, the 7 \q( 1410) meson, 
and with isospin 7 = 0, the /o-mesons at 980, 1370 
[1,2] and 1500 MeV from GAMS, CRYSTAL BAR¬ 
REL [1,2], WA102 [3] and BES [4]. Besides these reso¬ 
nances, there are different experimental indications [2], 
especially from BES [4], E791 [5], FOCUS [6], KLOE 
[7], SND [8], CMD2 [9], BELLE [10], WA102 [3] and 
7T7T scattering data [11,12] for some other scalar states, 
with 7 = 0, the er(600), / 0 (1710) and / 0 (1790), and 
with 7 = 1/2, the k( 840). The real quark and/or gluon 
contents of these states are not fully understood, which 
the interpretation using effective theories most of them 
based on a linear realization of chiral symmetry can¬ 
not clarify. In the following, we shall focus on the 
tests of the qq, qqqq and gluonium natures of these 
scalar mesons by confronting the recent experimental 
data with some QCD predictions based on QCD spec¬ 
tral sum rules (QSSR) complemented with some low- 
energy theorems (LET) [13-16] and lattice calculations 
[17,18], 

2. The 7 = 1, 1/2 scalar mesons 

The a 0 (980) and /«(840) masses 

These channels are expected to be simpler as we do 
not expect to have any mixing with a gluonium. If 
one assumes that these states are qq mesons, one can 
naturally associate them to the divergence of the vector 
currents: 

a 0 (980) -> = (m u - m d ) : u{i)d : , 

k(840) -> dpV? = (m u - m s ) : u(i)s : . ( 1 ) 

Within the QSSR approach, the properties and impli¬ 
cations of these mesons can be studied from the two- 


point correlator: 

= * J d 4 xe^(0|T^R/;(x)^V^(0)+|0). (2) 

Since the pioneering work in [19], numerous authors 
have used and improved the analysis of the previous 
correlator for the extraction of the running u-d mass- 
difference and strange quark masses [13,20]. The im¬ 
provements come from the inclusion of higher order 
terms in the PT QCD series [21,13]; the inclusion of 
the 1/q 2 term [22,23] which mimics the UV renormalon 
effects and which is also an alternative to the direct 
instanton effects where the later is not under a good 
quantitative control due to the incertainties of the in¬ 
stanton size and widths,...; the treatment of the spec¬ 
tral function using new Ktt phase shift data [24]. The 
obtained value of the u-d mass-difference is consistent 
with some other determinations, while the extraction of 
the ao(980) mass from the sum rule [13,20] is perfectly 
consistent with the data. Recent analysis [24] has lead 
to a value of the strange quark running mass which is 
consistent with some other QSSR determinations from 
different channels [25]. All these features support the 
qq meson assignements for the ao(980) and the ac( 840) 
mesons. However, the apparent wrong splitting of the 
ao(980) and the k( 840) mesons is intriguing. Here, we 
investigate this analysis using the ratio of exponential 
Laplace/Borel sum rules: 


Using the PT series to order o^, including the NP con¬ 
densates of dimension 6 and the new 1/q 2 terms, we 
give the prediction for the ao(980) in Fig. la) and for 
the k( 840) in Fig. lb) using a NWA. We use the most 
recent value of m s (2 GeV)=(96 ± 5) MeV compiled in 
[25]. One can see from these figures that the method 
reproduces the wrong splitting of the two mesons. The 
reason is that the SU(3) breaking effects increase the 


log / dt e * T Im ipud(t) . 


( 3 ) 


f 


1.1 1.2 1.3 1.4 



Figure 1. r in GeV 2 -dependence of the a) Mj( ao ) in GeV and b) M n in GeV for m s {2 GeV)=98 MeV at a given 
value of continuum threshold stability point: t c ~ 1 GeV 2 . 


value of the sum rule optimization scale compared to 
the one of the a o and then emphasize the contribution 
of the dimension six condensates, which is a vital cor¬ 
rection in the analysis. Indeed, from the analytic ex¬ 
pression of the sum rules, one can qualitatively extract 
the approximate mass formula 1 : 

M l - M l 0 + ~ 87t 2 to s (ss)to 

+ 2~sT pas ^ ^ ' T ° 

- ^M 2 r 2 r 0 , (4) 

where all different parameters including the ag mass 
are evaluated at the sum rule optimization scale r 0 — 
0.8GeV -2 ; p ~ 2 [29] indicates the deviation from 
the vacuum saturation of the four-quark condensate; 
( ss)/{uu) ~ 0.8 measures the SXJ{ 3) breaking of the 
quark condensate [13]. As shown in the above formula, 
the SU (3) breaking corrections are relatively small, and 
the four-quark condensate tends to decrease the k mass 
2 . Compared to the value of the n mass given in Fig. 
lb) using NWA, the finite width correction reduces the 
mass by about 20 MeV, where we have used the width 
of 310 MeV from BES [4]. One should mention that at 
the value of r 0 , the OPE in powers of r converges, while 
the radiative corrections to the parton graph, though 
large in individual sum rules, remain small in the ratio 
of moments 7 Z(r), as these corrections tend to com¬ 
pensate each others, justifying the uses of the result at 
this relativley low scale. From the previous analysis, 
we deduce: 

M ao ~ 930 MeV and M K ~ 920 MeV , (5) 

with about 10 % error, in good agreement with recent 
data [1,4,5] 3 . 

The decay constants 

The decay constant f ao of the clq normalized as: 

(01 W> 0 ) = ^faoK > ( 6 ) 

1 Notice that analogous formula in the vector channel explain 
with a relatively good accuracy the well-known (j)-p and K*-p 
mass splittings [26]. 

2 Using the instanton liquid model [27], the instanton contribu¬ 
tion has been explicitly shown in [28] to be negligible. 

3 Previous sum rule analysis of the k parameters [13] used as in 
put a value of m s and t c much higher than here and gives higher 
value of the k- mass. 


in the same way as f n = 92.4 MeV has been estimated 
several times in the literature [13,30]: 

f„ 0 — (1-6 ± 0.5) MeV , (7) 


where a better accuracy is claimed in [31]. Using SU (3) 
symmetry and the almost degeneracy of the oo and n 
masses, we expect to have with a good accuracy: 


f K _ rn s - m u 
fa 0 ~ m d - m u 


( 8 ) 


The hadronic couplings 

The ao and k hadronic couplings have been obtained 
using either a vertex sum rule [32,33] or/and 517(3) 
symmetry rotation [34]. The leading order vertex sum 
rule results are: 


9a 0 K+K- 


9a 0 K+K- 

9k,K+7T~ 


8ir 2 m s (ss) ( _ 
3x/2 M 2 K f K ^ 

e -( M K~ m l) T 0 (\ — 



~ 3 GeV 
-1.17, 


(9) 


where we have used [13]: m s (ss) ~ —0.8117^-/^, r = 
(ss)/(uu) — 0.8. and tq ~ 1 GeV -2 . We expect an 
accuracy of about 20% (typical for the 3-point fonc- 
tion sum rules) for these estimates. Using the 517(3) 
relation: 


9a 0 ryK — 


p 

Y 2 9a 0 K+K- 


one obtains: 


r (a 0 -> r?7r) 


\9oorjjp_ ( K 

167rM ao ^ Ml o 

84 MeV , 


( 10 ) 


( 11 ) 


in agreement with the range of data from 50 to 100 
MeV given by PDG [1]. Using the previous value of 
the k coupling, one can deduce: 

F(k Kn) ~ |r(« K+n-) ~ 104 MeV , (12) 

which is about a factor 4 smaller than the present data 
[1,4], but is a typical value for the width of a qq state. 

The 77 widths 

The 77 width of the oq has been evaluated using vertex 


















sum rules within the qq and four-quark assignements of 
this meson, with the result [30,34] : 


r a 0 (qq) - 

-7 77 ~ (0.3 ~ 2.1) keV , 

(13) 

and [30]: 

ra 0 (4g) - 

-> 77 ~ (2 ~ 5) x 10~ 4 keV , 

(14) 


to be compared with the data of (0.24±0.08) keV com¬ 
piled in [ 1 ], Due to the inaccuracy of the QSSR pre¬ 
dictions, we shall definitely use, in the following, the 
measured value of the do width as an input for our the¬ 
oretical predictions of the qq meson-width. 
Concluding remarks 

• The previous analysis shows that the mass and widths 
of the do (980) are well described by a qq assignement 
of this meson. 

• A QSSR analysis of the four-quark assignement [35] of 
these states gives predictions which reproduce the ex¬ 
perimental mass of the do(980) [36,30,13], like do the 
lattice calculations [37]. The result for the hadronic 
coupling g ao K+K- in the four-quark scenario depends 
crucially on the operators describing the do and can 
range from 1.6 GeV [38] to (5-8) GeV [30]. However, 
the prediction for the qn can agree with the data [30] 
depending on the size of the operator mixing param¬ 
eter. Therefore, an eventual selection of the two ap¬ 
proaches will be a precise measurement of g ao K+K~ 
or/and a lattice measurement of the decay constant 
which should depend linearly on the light quark mass 
in the qq scheme, but is a constant in the four-quark 
one. 

• Another problem arises when one computes the 77 
width using vertex sum rules. The ratio of the widths 
in the two approaches is [30] 4 : 

^4<m'/r<2977 ~ (1 ~ 2) X 10 , (15) 

which is of the order of ( a s /pi ) 2 , indicating that the 
four-quark assignement prediction is too small contrary 
to some claims in the literature. 

• For the n meson, the qq assignement can reproduce 
quite well its wrong splitting with the do, but fails to 
reproduce its large experimental width [about a factor 
4 smaller (see Eq. (12)], while the four-quark one gives 
a width of about a factor 2 smaller [38]. The failure 
of the two separate approaches (qq and four-quark as¬ 
signement) may suggest that the quark content of the 
k( 841) is more complex than naively expected: it can 
be a mixing between a qq and a four-quark states, or 
it can come from a large interference of the qq ground 
state with non-resonant backgrounds. Further tests are 
needed for clarifying its nature. 


3. The 1 = 0 bare scalar mesons 


The isoscalar scalar states are especially interesting in 
the framework of QCD since, in this anomalous U( l)y 
channel, their interpolating operator is the trace of the 
energy-momentum tensor: 


= jP( a s)G 2 + ^[1 + 1m(a s )]m i 'il>i'ipi , 

_ j 

4 We plan to come back to this point in a future work. 


where is the gluon field strengths, ^ is the quark 
field; (3(a s ) = (3i (a s / 7 r) + ... and 7 m (a s ) = 71 (a s / 7 r) + 
... are respectively the QCD /3-function and quark 
mass-anomalous dimension ((3\ = —1/2(11 — 2n/3) and 
71 = 2 for n flavours). In the chiral limit nii = 0, 9(( is 
dominated by its gluon component 0 g , like is the case of 
the 7 / for the U(1)a axial-anomaly, explaining why the 
?/-mass does not vanish like other Goldstone bosons 
for mi =0. In this sense, it is natural to expect that 
these 1 = 0 scalar states are glueballs/gluonia or have 
at least a strong glue component in their wave func¬ 
tion. This gluonic part of 9^ should be identified with 
the U(l)v term [39] in the expression of the effective 
lagrangian based on a U(3)l x U(3)r linear realization 
of chiral symmetry (see e.g. [40,41]). 

Unmixed 1 = 0 scalar qq mesons 

We shall be concerned with the mesons S 2 and S 3 

mesons associated respectively to the quark currents: 


J 2 = m : —=(uu + dd) : 
v 2 


and J 3 = m s : ss : 


(17) 


From the good realization of the SU( 2) flavour sym¬ 
metry (in u = nid and (uu) = (dd )), one expects a 
degeneracy between the a 0 and S 2 states: 

Ms 2 ~ M ao ~ 930 MeV , (18) 

while its hadronic coupling is [34,15]: 

9s 2 ir+n- ^ ^(uu)r 0 e^ ~ 2.46 GeV . (19) 

corresponding to 5 : 

r(S 2 -► 7 T+ 7 T-) ~ 120 MeV, (20) 

Using SU (3) symmetry, one can also deduce: 

gs 2 K+K- — ^gSiTT+ir- — 1-23 GeV . (21) 


The S 2 77 width can be deduced from the one of the 
ao(qq) obtained previously, through the non-relativistic 
relation (ratio of the square of quark charges): 

25 

Ts 2 ^77 - y r ao -77 - (0-7 ± 0.2) keV . (22) 

The mass of the mesons containing a strange quark is 
predicted to be [15]: 


M s JM k ~ 1.03 ± 0.02 => Ms 3 ~ 948 MeV , (23) 

if one uses M K = 920 MeV 6 , while its coupling to 
K+K~ is [15]: 


gs 3 K+K- - (2.7 ± 0.5) GeV . (24) 


These results suggest that the naive qq assignement of 
the er(600) = S 2 does not fit the data. 

5 We use the normalization: 


r(cr S 


■ 7T7r) 


3 \da B n+n- P A 4m^. \ ^ 
2 167 rM aB ^ 


(16) 


6 In [15] a higher value has been obtained because one has used 
as input the experimental mass Kq = 1430 MeV. 



Gluonia masses and decay constants 

• These states have been explicitly analyzed in [16,42, 
15] using QSSR of the two-point correlator: 

Ml 2 ) = 16 i J d 4 xe iqx ( 0 \T 0 ^(x) 6 £ (0)110) , (25) 

in the chiral limit m q = 0, for the observables: 

/»oo 

C n (r) = — dt t n e~ tT lmil; s (t) , (26) 

Jo 

and: 

7 ,00 

Rn,n+iW = l°g / dt t n e~ tT lm'i/j s (t) , (27) 

dr J o 

where n = —1, 0, 1, 2. For n = —1, the sum rule is 
sensitive to the subtraction constant: 

MO) = -16—(a s G 2 ) , (28) 

7r 

fixed from LET [43], where ( a s G 2 ) = (0.07 ± 0.01) 
GeV 4 [44,13]. One has found in [42] that, due to ^(O), 
the subtracted sum rule (SSR) £_i is more weighted 
by the high-energy behaviour of the spectral integral 
than the unsubtracted sum rule (USR) £o,i, 2 , which 
motivated the introduction of 2 resonances (below and 
above 1 GeV) for solving the controversial results ob¬ 
tained in the past. The results of the analysis using the 
standard OPE by retaining higher order PT series and 
the lowest dimension condensates are [16,42,15]: 

f(j B * (884 ± 116) MeV , 

fa - (390 ± 145) MeV , (29) 

corresponding to: 

M<j d ~ 1 GeV and M G ~ (1.5 ± 0.2) GeV. (30) 

M OB has been obtained in [16] using a least square fit 
of 7^o,i- Its decay constant has been obtained using a 
least square fit of IZq.i or/and a stability criterion of £o 
[42] 7 . The mass of the 2nd resonance has been fixed 
from 7?.o,l [15] and its decay constant comes from £_! 
[42]. We shall see that Mq is the one which can be com¬ 
pared with the present lattice value of about 1.6 GeV 
in the quenched approximation [17], while the erg mass 
will be a //-like meson expected to couple strongly with 
Goldstone boson pairs (huge OZI violation) [42,15] and 
playing a role in the saturation of the 17(1 )v two-point 
correlator subtraction constant ip s ( 0). It can only be 
tested using lattice with dynamical fermions (see e.g. 
[18] for the inclusion of the disconnected part of the 
scalar propagator). 

• One can also notice that a possible effect of the ra¬ 
dial excitation of the a can be obtained by matching 
the radial excitation contribution with the QCD con¬ 
tinuum. Assuming its mass to be around 1.4 GeV, one 
can deduce [42]: 

fa' B < (139 ~ 224) MeV , (31) 

7 A similar analysis in the U(l )a channel has given an estimate 
of the r\ parameters and of the U( 1) 4 topological charge [45,16]. 


while a weaker bound of about 500 MeV has been al¬ 
lowed in [15]. 

• The effect of the 1/q 2 term, to the previous results, 
which is an alternative of the direct instanton contribu¬ 
tion has been shown to be small [22], though this term 
is necessary for solving the sum rule scale hierarchy of 
the gluonia channels compared to the usual qq mesons. 

• A recent QSSR analysis of the same gluonium correla¬ 
tor using Gaussian sum rules and including instantons 
[46] confirms the previous mass values obtained in Eq. 
(30), but not the results in ref. [47,48], where it is ar¬ 
gued that the presence of the direct instantons solve 
the controversial results noticed in [42,15] between the 
subtracted n = — 1 sum rule with the other n > 0 un¬ 
subtracted ones , without the need of two resonances. 
In our normalization, the results in [47] are: 

M s = (1.25 ± 0.2) GeV, f s = (3 ± 0.3) GeV . (32) 

The mass value 8 does not contradict the ones in Eq. 
(30) as it is about the mean value of the two resonances 
ones, while the decay constant leads to [43,42,15]: 

-> 57 ) ~ 1.5 x 10“ 2 , (33) 

which is about 10 times higher than the one of the 
J/'0 —> /2(1.24)7 and which is already excluded by the 
BES [4] and some other data. This fact may signal 
some eventual internal inconsistencies in the treatment 
of the instanton contributions. 

• An upper bound on the gluonium mass has been also 
obtained in [15]: 

M g < (2.16 ± 0.22) GeV , (34) 

using the positivity of the moment 7^1,2- This bound 
has been strengthened in [47] to 1.7 GeV, which cannot 
be an absolute bound because of the inclusion of the 
QCD continuum model and of its related uncertainties 
for its derivation. 

Gluonia widths to 7T7 t 

• For this purpose, we consider the gluonium-7T7 t vertex: 

V(q 2 ) = (n 1 |0£|7i-2), q = Pi~P2, (35) 

where: V(0) = 2 m 2 . In the chiral limit ( m 2 ~ 0), the 
vertex obeys the dispersion relation: 

F( " 2) = l 7 ImF(() • (36) 

which gives the 1st NV sum rule [42]: 

\ ]T g s ^V2f s ~0. (37) 

S = CTB 

Using the fact that V'(0) = 1 [49], one obtains the 2nd 
NV sum rule [42]: 

4 gSmrVZfs/Mg = 1 . (38) 

S=ub ,<y' B ,G 

To a first approximation, we assume a crg-dominance in 
the 2nd NV sum rule, while in the 1st sum rule, there is 

8 It is an upper bound in [48]. 



a matching between the ap and all higher mass glueball 
contributions, which we replace by an effective o' mass 
of about 1.4 GeV. Then, one obtains: 

9a B TT+TT- — 9a B K+K- — (3.2 ~ 6 . 8 ) GeV . (39) 

A complete matching in the 1st NV sum rule would 
lead to an effective coupling: 

9a> B *+, r— « 12 GeV , (40) 

when using f' a ~ 500 MeV. An interpretation of 
this value of g a ' B ^ is unclear but it is expected to 
parametrize all higher states contributions to the 1 st 
sum rule. If one uses a resonance mass of about 1.4 
GeV, one would obtain a very broad a' B which can 
mimic the red dragon proposed earlier [50], but can¬ 
not manifest as a peak, making its identification with a 
true resonance difficult. Different strategies for extract¬ 
ing the /o(1370) parameters from the continuum back¬ 
ground has been discussed in [ 12 ] leading to a width of 
the fo compiled in PDG [1]. For definiteness, we shall 
use in the following, the value of —> 7 T 7 t about 

250 MeV obtained in [12] within a factor two accuracy. 
This leads to: 

9^*+*- - 9a' B K + K- ^ (4.2 ± 0.7) GeV . (41) 

ctb as an ^Mike meson 

• One can check that the LET for if> 3 (0) given in Eq. 
(28) is almost saturated by the contribution of the low¬ 
est mass <tb{ 1 .) meson: 

V’s(O) — 2 Mj B /J B ~ 1.5 GeV 4 , (42) 

compared to the LET value of 1.6 GeV 4 . This prop¬ 
erty is very similar to the contribution of the rf in the 
topological charge (subtraction constant of the anoma¬ 
lous U(1)a gluonium two-point correlator [51]),which 
explains why it is not degenerated with the pion at fi¬ 
nite N c . 

• The as large coupling to pseudoscalar pairs, through 
OZI violating process, can also be compared with the 
affinity of the rf to couple to ordinary mesons, making 
it as an ambidextre gluonium-meson state 

• Like in the (7(1) a sector, the quenched lattice simu¬ 
lations obtain a higher glueball mass which is not the 
rf mass. The rf - mass and decay constant are only mea¬ 
sured from a lattice calculation of the (7(1) a topolog¬ 
ical charge [52], We expect that the same situation 
occurs in the U(l)v channel, where a quenched lattice 
gives a scalar gluonium mass of about 1.5 GeV, while 
the <jb parameters can be obtained from the measure¬ 
ment of the scalar correlator subtraction constant ip a ( 0 ) 
including dynamical fermions. 

G(1.5) widths into r/'r/', r]ri' and r/r] 

• Analogous low-energy theorem [42] gives: 

(mK\m) = 2M^, (43) 

where 771 is the unmixed 17(1) singlet state of mass 
M m ~ 0.76 GeV [51]. Writing the dispersion relation 
for the vertex, one obtains the NV sum rule: 

4 E 9s mm V2f s = 2 M^, (44) 

S=a B ,a' B ,G 


which, by assuming a G-dominance of the vertex sum 
rule, leads to: 

g Gmm » (1.2 - 1.7) GeV. (45) 

Introducing the “physical” rf and 77 through: 
rf ~ cosOprji — sin 6 pr]$ 

77 ~ sin 0 p 77 i + cosOprjs, (46) 

where [1,53] Op ~ —(18±2)° is the pseudoscalar mixing 
angle, one can deduce: 

r(G -7 rfrf) ~ (5 ~ 10) MeV. (47) 

The previous scheme is also known to predict (see NV 
and [59]): 

r = - 0-22, 9g vv - sin 0 P g G r,r,', (48) 

J- Grjry 

compared with the GAMS data [1] r ~ 0.34 ± 0.13. 
This result can then suggest that the G(1.6) seen by 
the GAMS group is a pure gluonium, which is not the 
case of the particle seen by Crystal Barrel [1] which 
corresponds to r ss 1 . 

Gluonia widths into 4-7T 

Within our scheme, we expect that the 47 r are mainly 
S'-waves initiated from the decay of pairs of <tb. Using: 

{vbW\<tb) = 2MZ b , (49) 

and writing the dispersion relation for the vertex, one 
obtains the sum rule: 

\ E 9Sa B , B V2f s = 2Ml B . (50) 

i—a B ,u' B ,G 

We use M„ b ~ 0.6 ~ 1 GeV, M a < s ~ 1.4 GeV and the 
observed /o(1.37) width into 47 t of about (106 ~ 250) 
MeV [1] (S-wave part). Neglecting, to a first approx¬ 
imation, the <jb contribution to the sum rule, we can 
deduce: 

\9Ga B cr B \ ~ 1-3 ~ 3.7 GeV , (51) 

where the first (resp. second) value corresponds to 
M tJB ~ 0.6 GeV (resp- 1 GeV). This leads to the width 
into cr(600)cr(600) of about (7 — 55) MeV, much larger 
than the one into 7777 and 7777 k This feature is satisfied 
by the G(1.5) state seen by GAMS, Crystal Barrel and 
WA102 [1], However, the previous approaches show the 
consistency in interpreting the G(1.5) seen at GAMS as 
an “almost” pure gluonium state (ratio of the r\rf ver¬ 
sus the 7777 widths), while the state seen by the Crystal 
Barrel and WA102, though having a gluon component 
in its wave function, cannot be a pure gluonium be¬ 
cause of its prominent decays into 7777 and 7 r + 7 r~. 
Gluonia widths into 77 

These widths have been derived in [42] by identifying 
the 77 -glue-glue box diagram with the scalar 77 La- 
grangian where the quarks in the internal have been 
taken to be non-relativistic. In this way, one has ob¬ 
tained: 

5ct77 cs ^V 2 f a Ml E ’ 

^ ' u.d.s 


( 52 ) 



where Qi and ?Tij are the quark charge and constituent 
masses. This leads to: 

T (a B -> 77 ) ^ (0.03 ~ 0.08) keV , 

T (a' B -> 77 ) ~ (0.01 ~ 0.03) keV , 

T(G -► 77 ) ~ (0.3 ~ 0.6) keV . (53) 

Alternatively, one can use the trace anomaly to order 
k 4 in order to deduce: 

< 0 |i/ 3 MG 2 | 7 i 72 > ^ -<0| (^) RFrFr\I 1 I 2 ) ,(54) 

where R = 3J±Q?, ct is the QED coupling. This 
relation gives: 

fairy - (^) R • ( 55 ) 

(7,... 

From this relation and using the previous values of f a 
and f' a , one can deduce: 

T G ^ 77 ~ (1 ~ 6 ) keV , (56) 

which is quite inaccurate but still consistent with the 
previous determination. 

Comments 

Comparing the above results, especially the predicted 
widths, with the experimentally observed candidates, 
it is likely that the er and some of its radial excitations 
have a lot of glue in their wave functions. As a conse¬ 
quence, a quarkonium-gluonium decay mixing scheme 9 
has been proposed in the 1 = 0 scalar sector. [34,15], 
for explaining the observed spectrum and widths of the 
possibly wide er(< 1 GeV) and the narrow /o(980) 10 . 

4. Meson-gluonium mixing below 1 GeV 

• BES data suggest that the cr(600) is produced in the 
OZI forbidden J/tp —> (pn + ir~ process [4], which can 
indicate the large amount of glue in its wave func¬ 
tion. Its production from the OZI allowed J/tp —> 
W 7 r + 7 r _ , K + K~ processes, is expected to be due to 
its quark component, while the relative small branch¬ 
ing ratio in the OZI allowed J/tp —> (pK + K~ process 
relative to J/ip —■» (pir + ir~ can be due to an interfer¬ 
ence between the K + K~ amplitude from the gluon and 
quark components of the a. 

• In the same way, the /o(980) is produced in the OZI 
violating J/xp —> <pir + n ~, [4] and D s —> 7 r~ 7 r + 7 r + [5] 
processes which may also indicate its gluonium compo¬ 
nent, while its production from J/vp —> (pK + K ~, can 
signal a strong ss component in its wave function. We 
shall keep in mind these results for building the mixing 
scheme. 

The meson-gluonium mixing scheme 

• We assume that the observed states come from the 

9 This has to be contrasted with the small mass-mixing coming 
from the off-diagonal two-point function [54]. 

19 We shall not consider in our analysis the recent result of [55] 
from 77 —> AW where the resulting a —» 77 width of about 4 
keV is much bigger than generally expected and which needs to 
be confirmed by some other data. We plan to come back to this 
point in the future. 


mixing between the gluonium as and quark S 2 = 
\/\/ 2 {uu ± dd) and S 3 = ss bare states 11 : 

a \ _ f cos 0 s sin 6 s \ ( cr B \ 
fo )~\ — sin 9 S cost's ) V S 2 + <PsS 3 ) 

where (ps = 1 /x /2 for a SXJ{ 3) singlet and —y/2 for an 
SU(3) octet. 

• In [34], the mixing angle 0s has been fixed from 
the analysis of the predicted decays of the hypothet¬ 
ical bare states S?, a B and of the observed meson 
/o(980) —> 77 . Using the predictions: r(5 2 —> 77 ) 
in Eq. (22), T (a B —> 77 ) ~ 0.03 keV [42], and the 
recent data T(/ 0 (980) —> 77 ) ~ (0.4 ±0.1) keV [1], one 
can deduce: 

9s ^ (45 ± 15)° , (57) 

where the S 3 —> 77 width is suppressed as (2/25) com¬ 
pared to the one of the 5 2 . 

• Using, for definiteness, as inputs the theoretical pre¬ 
dictions given in Eqs (19) and (39), we predict the cou¬ 
plings: 

9**+*- - (5.3 ±1.8) GeV , 

9fo*+*~ - (1.8 ±1.3) GeV, (58) 

in reasonnable agreement with the ones from the data: 

CV - (3.0 ±1.5) GeV, 

9 e fZ+«- - (1-5 ±0.3) GeV, (59) 

corresponding to: 

r exp (c7(600) 7 T 7 r) « 481 MeV , 

r e ^(/ o (980) 7 T 7 r) w (70 ± 30) MeV , (60) 

• In order to predict the mixing parameter <pg, we fit 
the experimental foK + K~ width from BES [4]: 

9f 0 K+K-/9f 0 -K+n- - 2-05 ± 0.15 , (61) 

and we use the theoretical predictions given in Eqs (19), 
(24) and (39). Then, we obtain: 

(ps — 3.0 , and 9 &K+ k- / 9 ctk+ 7 t- — 2 . (52) 

Further data are needed for improving and testing this 
result. 

• One should note that using only the constraint from 
the fo —* 7 r + 7 W and fo —> K + 1\~ hadronic widths, one 
would obtain: 

0s = 16° , and (ps — —1-4 , (63) 

indicating that the a is an almost pure gluonium and 
the fo a qq SU( 3) octet, which is similar to the scheme 
in [50]. However, the result in Eq. (63) would give a 
too high value of fo —■► 77 and does not explain the OZI 
violating production of the fo in J/ip and <p radiative 

11 In our approach, we first calculate the real part of the masses 
of these hypothetical states and deduce their widths using vertex 
sum rules. These states would correspond to the bare states in 
the A'-matrix formalism (see e.g. [12]). Due to the large error in 
our mass predictions, we neglect, to a first approximation, some 
possible shifts on the masses which can be induced by the decay 
processes mentioned in [12]. 



decays. 

Comments on alternative approaches 

• One should note that a four-quark QSSR analysis 
gives [38]: 

9f 0 K+K- — 9a 0 K+K- — (1-6 ±0.1) GeV 

<7/0*+*- ~ (0.47 ± 0.05) GeV , (64) 

where the absolute values differ from the ones given 
in [56] and then question the realibility of the results 
obtained there. Therefore, lattice calculations of these 
couplings become mandatory. Eq. (64) leads to: 

r(/ o (980) -4 tttt) « 7 MeV , (65) 

which is too small compared with the range (40 ~ 100) 
MeV given by the data [1]. 

• Alternative approaches based on KK loop, KK 
molecules and four-quark states can predict value of 
F(/o(980) —> 77 ) in agreement with the data [57],[58], 
but, most of them, do not give a satisfactory predic¬ 
tion for the fo —► 7r + 7r - width and their production 
from OZI violating decays. On the other, it would be 
interesting to see the connection of these effective ap¬ 
proaches with the quark-gluon picture used here. 
Tests from J/ip and <fi radiative decays. 

• These decays are known to be a gluonium filter. The 
production of a gluonium S from J/ip radiative decays 
can be approximated by [49]: 

r(j/v.- 7 S)= !s gL_(^f) 1 (|!f) 4 

fs > ( 66 ) 

where M c ~ 1.5 GeV is the charm constituent quark 
mass. In our scheme, the a is mostly a gluonium. 
Therefore, one expects the branching ratio: 

B{J/ip -4 7 cr) « 19 x 1CT 5 . (67) 


(l-Ml/M^)' 
r >e+e-) 


an isoscalar are CKM-suppressed due to the c-u tran¬ 
sition at the weak vertex. 

• Scalar gluonium or/and ss productions 

Semi-qualitative but model independent results for the 
production of gluonium have been given in [28] (see 
also [61]): 

— The only way to obtain a non-CKM suppressed 
isoscalar is to look at the semileptonic decay of the 
ZU-meson, where the light quark is a strange one and 
an isoscalar ss or/and gluonium state can be formed. 

— If the ss state is relatively light (< 1 GeV), which 
might be the natural partner of the (uu + dd) often 
interpreted to be a a (600) in the literature, then, one 
should produce a KK pair through the isoscalar ss 
state. The non-observation of this process will dis¬ 
favour the qq interpretation of the cr and fo mesons. 

— If a gluonium state is formed it will decay with even 
strength into 7T7 t and a KK pairs. Therefore a gluo¬ 
nium formation in semileptonic D s decays should result 
in the decay patterns: 

D s — » cref± — > 7T7rf± D s — > <jb(-v —► KKtu , (71) 


with about the same rate up to phase space factors. 
The observation of the semileptonic 7T7 t decay of the 
D s by E791 [5] is a sign for glueball formation 12 . 

— Using, e.g., the result in [60], the one for light S(qq ) 
quarkonium production behaves as: 

T[D S -4 lu] ~ |V C9 | 2 G^M c 5 |/ + (0)| 2 . (72) 


— The <t b (gg) production, can be obtained from the 
1 /M c behaviour of the WWgg box diagram. Using 
dispersion techniques similar to the one used for J/ip 
a B J processses [43,42,13], one obtains, assuming a D s 
and a B -dominances [28]: 


T[D S -4 <j B lu] 


\V cs \ 2 G 2 f 


m<PsG 2 \cr B )\ 2 

M C M± 


(73) 


Extrapolating the previous exppression to the (p -meson 
and using M s ~ 500 MeV, one obtains: 


where (0|</>sG 2 |ctb) is by definition proportional to 
/ ct M 2 . Then, one deduces: 


B{cp —> 7 <t) « 12 x 10 -5 , (68) 

which, despite the crude approximation used, compares 
quite well with the KLOE [7] data: 

B(cp -4 77r°7r°) « (10.9 ± 0.3 ± 0.5) x 10“ 5 . (69) 

Tests from semileptonic decays 

This section has been discussed in details in [28] and 
will be only sketched in the following. 

• S 2 (uu + dd) meson productions 
If the scalar mesons were simple qq states, the semilep¬ 
tonic decay width could be calculated quite reliably us¬ 
ing QSSR, where the relevant diagram is a quark loop 
triangle. Several groups [60] predict all form factor to 
be: /+(0) « 0.5 , yielding, for Ms 2 — 600 MeV, a decay 
rate: 

T(D -4 S 2 lv) = (8 ± 3)10 -16 GeV , (70) 

which is, unfortunately, even in high stastistics experi¬ 
ments, at the edge of observation since the decays into 


T{Ds -> a B {gg ) lv\ 1 ( fa _\ 2 (74) 

T[D S —> S q (qq) lu] ~ |/+(0)| 2 \M C ) ’ { ’ 

which is 0(1) for f„ ~ 0.8 GeV. 

• Comments 

These semi-quantitative results indicate that the gluo¬ 
nium production rate can be of the same order as the qq 
one contrary to the naive perturbative expectation (a 2 
suppression), which is a consequence of the OZI-rule 
violation of the a B decay. 13 However, it also shows 
that, due to the (almost) universal coupling of the <7 B 
to Goldstone boson pairs, one also expects a produc¬ 
tion of the KK pairs, which can compete with the one 
from ss quarkonium state, and again renders more dif¬ 
ficult the identification of the such ss state if allowed 
by phase space. 

12 An alternative explanation assuming fo(ss) and using KK 
loop has been given in [58]. However, the same assumption for 
fo but using QSSR leads to a negative conclusion [62], 

13 Productions of the scalar mesons in B-decays have been dis¬ 
cussed in [63]. 



5. Properties of the mesons above 1 GeV 


where P = n, K ; n ~ 2. Then, we deduce: 


QSSR does not have a precise systematic framework 
for extracting the properties of the radial excitations, 
except the approximate value of the mass indicated by 
the value of the QCD continuum threshold at which the 
mass of the ground state has been optimized. In order 
to check this result, one often uses (if available) empiri¬ 
cal observations of the splittings between the radial ex¬ 
citations and the ground state, or/and arguments based 
on the (linear) Regge trajectories. For this reason, the 
discussions which we shall give below will be very qual¬ 
itative. 

The 1=1 and 1/2 mesons 

The ao(1450) and the K/ (1430) are almost degenerated 
indicating the restoration of the SU( 3) flavour sym¬ 
metry where the SU( 3) breakings behaves like m 2 s /M B 
and m s (ss)/Mft for the radial excitations. These scalar 
radial excitations are also almost degenerated with 
the pseudoscalar 7t(1300) and K( 1460) also indicat¬ 
ing the restoration of the spontaneous breaking of the 
SU (3) l x SU ( 3)u chiral symmetry broken by the quark 
condensate at higher scale 14 . Within such observa¬ 
tions, one can expect that these scalar states are qq 
states 15 . 

The I = 0 scalar mesons 

There are proliferations of these states from the data 
[1,4]: / o (1370), / 0 (1500), / o (1710) and / o (1790). Us¬ 
ing the previous symmetry restorations, one may ex¬ 
pect that the 1st radial excitation of the (uu + dd) 
and ss unmixed states should be in the range of 1400 
MeV, which are in the vicinity of the /o(1370) and 
/o(1500) being good experimental candidates. The 2nd 
radial excitations are expected to be in the range of 
the 7r(1800) where good experimental candidates are 
the /o(1710) and / o (1790). Within these observations 
in mind, we shall interpret the different data given by 
BES [4]: 

• The /o(1370) and /o(1710) are respectively produced 
through the OZI forbidden J/ip —> pn + i r~ and al¬ 
lowed J/ip — * u>K + K~, processes [4], while the lat¬ 
ter is also produced from the glueball filter process 
J/ip —> r )K + K^. These features can indicate that they 
can have an important gluonium component in their 
wave function. To a first approximation and iteration, 

we assume that these two states result from the mixing 
16 . 


( /o(1370) \ ( cos 0' s sin 9' s \ ( a' B \ 

'v /o(1710) ) { -sin 0' s cosO's J \ ^ + ^ J 

which is a replica of the mixing among ground states 
in Eq. (57). We fix the coupling of the qq radial ex¬ 
citations by assuming that they behave like the one of 
the pion and 7r(1300), namely: 


9s^ 3 pp 


( M sL,X 

\mX) 9s2 ’ 3Pp 


(75) 


14 Similar though qualitative arguments have been given in [64]. 

15 Similar conclusions have been also independently reached in 
[65] using FESR. 

16 Some other alternative mixing schemes above 1 GeV have been 
also proposed [12,15], [66,50], [67,68]. 


9s' 2 ir+ir- — 2gS' 2 K+K- — 4.8 GeV , 

9S'K+K- - 6.4 GeV. (76) 

Using the experimental input from BES [4]: 


F(/o(1370 - 7T7T) 
r(/o(1370 -> KK) 
r(/o(1370 -► 7T7T) 


~ (265 ± 40) MeV, 
~ (0.08 ±0.08), 


(77) 


one can hx the two angles to be: 

e’ s ~ 75° <p’ s ~ -0.55 . (78) 

The solution indicates that the /o(1370) contains more 
(uu + dd) than a gluonium, while the /o(1710) contains 
more glue and ss than (uu ± dd), which is in line with 
general trends. From this analysis, we predict: 

9 f 0 (i7io)7T+7T- ~ 2.8 GeV , 

9fo(i7io)K+K- ~ 4-3 GeV , (79) 

which lead to the widths: 


r[/,f (1710) ->tt+tT ] w 90 MeV , 
r[/<f (1710) -► K+K~] « 176 MeV . (80) 

The results are of the order of the BES data, which are 
respectively < 16 MeV and 125 MeV [4], but are not 
yet quite satisfactory. As the /g ?l (1710) contains more 
glue after the mixing, it is natural that it will mix with 
the glueball G(1500) with the theoretical parameters 
obtained previously. 

• Therefore, in the 2nd step analysis, we consider that 
the observed / o (1710) and / o (1500) result from the 
mixing between f$ l ( 1710) and G(1500). 

• The /o(1500) is produced in the gluonium filter pro¬ 
cess J/ip —> 77r + 7r”, and it is needed for improving 
the fit of the OZI violating J/ip —> pm + Tr~ process, 
but has a small effect in the J/ip —> ifiKK process [4]. 
On the other, one expects from previous sections that 
the G(1500) couplings to 7T7 t and KK are small as the 
vertex sum rule is almost saturated by the <jb, but we 
do not have a precise quantitative control of such cou¬ 
plings. However, one can check that the model cannot 
explain simultaneously the small /o(1500) —» K + K~ 
and /q( 1710) —> 7r + 7r~ widths. 

• A solution to this problem may be given by the 
mixing of the previous two states with the /o(1790) 
and some 2nd radial excitations of the qq states ex¬ 
pected to be in this range of energy 17 . The /o(1790) 
may contain more gluon in its wave function as it is 
produced in the glueball filter channel J /ip —> 77r + 7r~ 
and in the OZI violating J/ip —> (/>7T + 7r~ process, 
while it also decays into 47 t. The smallness of the 
/o(1790) —» K + K~ width may come from a destruc- 

17 Alternative explanation has been in [69], where it is argued 
that the glueball /o(1710) has chiral coupling to pairs of Gold- 
stone bosons using perturbative QCD arguments. However, the 
validity of this result has been questioned in [70] due to instanton 
contributions. 



tive interference between the gluonium and qq states 18 . 


6 . Conclusions 

Motivated by the various recent experimental pro¬ 
gresses for producing scalar mesons, we have pro¬ 
posed new results and updated previous predictions of 
[15,42,34]: 

• The unusual wrong splitting between the ao(980) and 
k( 900) being respectively a ud and us mesons can be 
understood from QSSR using the present value of rh s . 

• The maximal meson-gluonium mixing for the cr(600) 
and /o(980) appears to be supported by the present 
data. 

• The ao(ud), k(us), <t( 600) and /o(980) mesons ap¬ 
pear to complete the scalar nonet below 1 GeV. 

• The productions of the 1=0 mesons above 1 GeV 
through the OZI violating J/ij) — <j)ir + ir~, D s —> 37t 
decays and J/tj) —> jS glueball filter processes may 
indicate that most of these 1 = 0 mesons have impor¬ 
tant gluonium in their wave functions, where a simple 
meson-gluonium mixing scheme can explain the gen¬ 
eral features of the data. 

• Our results suggest that the /o(1370) is mostly qq , 
while the /o(1500), /o(1710) and /o(1790) have sig¬ 
nificant gluonium component in their wave functions. 
Tests of these results can be provided by the measure¬ 
ments of the pure gluonium rj'r/ and 47 t specific U(1)a 
decay channels. 

We conclude that, after about a 1/4 century study, we 
still remain with more questions than answers on the 
true nature of scalar mesons. 
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